BI-TP 2012/49 
HIP-2012-28/TH 



J2' 

9-; 
en 



Lattice study of an electroweak phase transition 

at m h ~126 GeV 

M. Laine a , G. Nardini b , K. Rummukainen c 

^ Institute for Theoretical Physics, Albert Einstein Center, University of Bern, 
Sidlerstrasse 5, CH-3012 Bern, Switzerland 

h 'Faculty of Physics, University of Bielefeld, D-33501 Bielefeld, Germany 

c Department of Physics and Helsinki Institute of Physics, 
University of Helsinki, P.O.Box 64, FIN-00014 Helsinki, Finland 

Abstract 

We carry out lattice simulations of a cosmological electroweak phase transition for a Higgs 
mass m/i ~ 126 GeV. The analysis is based on a dimensionally reduced effective theory for 
an MSSM-like scenario including a relatively light coloured SU(2)-singlet scalar, referred to 
as a right-handed stop. The non-perturbative transition is stronger than in 2-loop pertur- 
bation theory, and may offer a window for electroweak baryogenesis. The main remaining 
uncertainties concern the physical value of the right-handed stop mass which according to 
our analysis could be as high as ~ 155 GeV; a more precise effective theory derivation 
and vacuum renormalization than available at present are needed for confirming this value. 
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1. Introduction 



One of the most important cosmological conundra is that we live in a Baryon Asymmetric 
Universe (BAU), meaning that very few antiprotons are seen in cosmic rays, although quarks 
and antiquaries are produced in almost equal measure at collider experiments. One of the 
possible explanations for BAU is called electroweak baryogenesis pQ: it makes use of the 
Higgs mechanism which is assumed to undergo a "first order" phase transition in the Early 
Universe. A phase transition of this type leads to thermal non-equilibrium, one of the nec- 
essary Sakharov conditions for explaining the BAU. The others (C, CP, and baryon number 
violation) are also part of the Standard Model and its simple extensions. The open issue, for 
any model, is whether these necessary conditions also amount to sufficient ones. 

Although electroweak baryogenesis is by no means the only available scenario for generating 
the BAU, it is attractive in that it offers for a very restricted framework, being successful only 
for specific parameter values that can conceivably be probed at the LHC. A disadvantage is 
that a complete and reliable calculation of the BAU is theoretically demanding (for a recent 
review see ref. [2]). However, a crucial step, which by itself leads to strong constraints, is 
to prove the existence of a strong first-order electroweak phase transition. For this step, 
which is the topic of the present paper, uncertainties can be brought under control through 
a combination of analytic computations and large-scale lattice Monte Carlo simulations. 

Lattice studies of variants of the electroweak phase transitions have a long history by now. 
They represent the only known systematic way of circumventing the infrared problem of 
thermal field theory [3] that limits the accuracy of perturbative evaluations. For instance, it 
was originally envisaged that electroweak baryogenesis could work even within the Standard 
Model pQ, but nowadays it is known that this possibility is not realized in nature. An un- 
ambiguous reason is furnished by lattice simulations [U [5] : unlike suggested by perturbation 
theory, the transition is a crossover for a Higgs mass compatible with either LEP or LHC [6l [7j 
bounds, which implies that the system does not deviate from equilibrium 

On the other hand, simple extensions of the Standard Model, particularly the Minimal 
Supersymmetric Standard Model (MSSM), might change the picture. Indeed, even though 
very strongly constrained by now [£l [10], realizing electroweak baryogenesis in the classic 
MSSM |llj-[20j (or in extensions of the Standard Model resembling MSSM at low ener- 
gies [21]) seems not to be excluded [22j [231 121] • In particular, if right-handed stops are 
sufficiently light and left-handed stops are heavy enough, the electroweak phase transition 
can be strong even for a (lighter CP-even) Higgs mass ~ 126 GeV [221 125] . This assertion 
relies, however, on perturbation theory, whose accuracy cannot be taken for granted. There- 
at is commonly believed that the Standard Model also contains too little CP violation to allow for elec- 
troweak baryogenesis, however this issue is harder to prove beyond doubt, cf. e.g. ref. [8] for recent work and 
references. 
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fore, to confirm or rule out the scenario, lattice Monte Carlo simulations appear welcome. 

In the past, lattice analyses have extensively studied the electroweak phase transition within 
the MSSM. However, even the latest simulations [26] only focussed on a part of the MSSM 
parameter space that corresponds to a lightest Higgs mass < 115 GeV. This seems to be 
excluded by the recent LHC data [27]. It thus appears well-motivated to repeat the lattice 
analysis for ~ 126 GeV, and this is the aim of the present study. 

The paper is organized as follows. In sec. [2] we review the dimensionally reduced effective 
theory of the MSSM with a light right-handed stop, and choose a parameter region where 
perturbative calculations find a strong electroweak phase transition. Within this region, we 
select a parameter point with ~ 126 GeV where we analyze the phase transition on the 
lattice (sec. [3]). The lattice results are compared with the perturbative ones in sec. HI Finally, 
sec. [5] is dedicated to an outlook and conclusions. 

2. Dimensionally reduced effective theory 

For the benefit of an impatient reader, we start by summarizing the four-dimensional (4d) pa- 
rameter values which the effective theory simulations are believed to correspond to (sec. 12. 1| ). 
Subsequently the theoretical foundations and practical implementation of the dimensionally 
reduced (3d) effective theory construction are briefly reviewed (sec. 12. 2\i . 

2.1. Parameter values 

The present study is based on lattice simulations as explained in ref. [26], and on analytic 
dimensional reduction and vacuum renormalization formulae as described in ref. [28] , The 
physics setting is different from what might appear ideal from today's perspective: in par- 
ticular, the presence of a relatively light CP-odd Higgs mass m A w 150 GeV, as well as 
of relatively light gluinos of mass M 3 < 300 GeV, were assumed in ref. [28] . Even though 
this setting is problematic because of light stop bounds (forbidding small M 3 [29]) and dark 
matter constraints (disfavouring light m A [30]), it provides for a conservative framework to 
prove the existence of a strong first-order phase transition with ~ 126 GeV. Indeed, we 
expect to find a stronger transition for larger CP-odd Higgs mass (decreasing m A weakens the 
transition [26} l3l] ) . Having relativistic gluinos in the thermal bath increases the right-handed 
stop thermal mass by ~ 20% and consequently pushes the right-handed stop vacuum mass 
parameter to more negative values but, as we have checked by a resummed 1-loop estimate, 
does otherwise not significantly affect the phase diagram. 

A major difference with respect to ref. [28] where a left-handed stop mass m Q < 1 TeV 
was assumed, is that here we push m Q to much larger values. This is needed to achieve 
rrih ~ 126 GeV but, as pointed out in ref. [32], this induces large logarithms that were not 
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resummed in ref. [28J. Our evaluations of the Higgs and right-handed stop masses, and 
m ( - fl , as functions of the MSSM parameters, are therefore approximate. 

Another important point is the renormalization scale at which the couplings appearing in 
the thermal mass corrections are evaluated. Following ref. [17], it was assumed in ref. |28| 
that the couplings run to a scale ~ 2irT like in the Standard Model [33] . To crosscheck 
the argument requires carrying out a full 2-loop dimensional reduction computation, which 
is absent at present in the parameter range considered. The concrete consequence of the 
assumption of a scale 2ttT ^> m top is that the strong gauge coupling is quite small; this 
implies that the right-handed stop squared mass parameter, mfj, does not need to be as 
negative as sometimes assumed; and subsequently, that the physical right-handed stop mass 
mi R ~ yj m^ op — mfj (for At ~ and fh% = —m 2 v > 0) is larger. 

To be reminded of these uncertainties, we tag the parameters mentioned, as well as the 
temperature, by a star in the following: 

m a , m Q , m h , m iR , T , (2.1) 

and similarly for the less significant parameters. These numbers are therefore not to be 
interpreted as precise physical values. 

With these reservations, the dimensional reduction formulae and the notation are identical 
to refs. \2Q\ 128], (As a small point, explicit CP violation has been switched off for simplicity). 
The perturbative phase diagram is illustrated in fig. (TJ For the parameter setting 

m* = 70.5 GeV , m* = 7 TeV , // = M 2 * = m* = 150 GeV , 

tan /3* = 15 , A* t /m* Q = 0.02 , (2.2) 

where //*, M|, tan/3* and A* are defined in a standard way (see e.g. ref. [28]), the pertur- 
bative calculation yields v(T*)/T* = 0.9 in Landau gauge, where v(T*) is the gauge-fixed 
expectation value of the lighter Higgs at the critical temperature T* (the precise definitions 
of these observables are given in sees. [3] and [4]) §| Notably, at this parameter point the lightest 
Higgs pole mass is m* h ~ 126 GeV within the 1-loop approximation [34] . 

Comparing with refs. |22} 125]. the mass parameter fh* u and critical temperature T* evalu- 
ated at nig = 7 TeV are substantially smaller, and consequently the physical stop mass m-^ 
is larger. The difference seems to be related to running effects in the thermal mass correc- 
tions as mentioned above (the discrepancy is smaller at m* Q = 1 TeV, but of course then m* h 
becomes unphysically light). It is our ultimate goal to improve the dimensional reduction 
and vacuum renormalization computations for the case of a very large nig, but unfortunately 
this requires a substantial amount of new work, which is postponed to future. If it turns 
out that our tagged parameters are close to the physical ones, in particular mr ^155 GeV 

2 The vev normalization is v(0) ~ 246 GeV. 
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Figure 1: Slices of the pcrturbativc phase diagram of the effective theory considered, in terms of the 
MSSM-likc parameters defined in the text. The upper boundary corresponds to parameters at which 
v(T*)/T* equals 0.9 in the Landau gauge; the lower boundary to values below which the theory is 
driven to a colour-broken minimum according to the 2-loop effective potential. 



(cf. fig- [Tj) , then the tension between LHC data and the considered MSSM scenario would be 
relaxed [9j \TU[ [22| [29] . (Note that since the transition we find is comfortably strong enough 
for baryogenesis, there remains some tolerance for a minor error in m i R -) 

Despite the above uncertainties, there is still a well-defined problem to solve. Namely, we 
compare non-perturbative Monte Carlo simulations with 2-loop perturbation theory within 
the 3d effective theory. In other words, the effective parameters entering the simulations and 
the analytic formulae are identical. Therefore the finding of a non-perturbative effect in one 
or the other direction is likely to also persist with slightly modified parameter values§| The 
actual simulations are carried out for the parameter values of eq. ()2.2|) . 



2.2. Review of theoretical setting 

In order to appreciate the methods to be used, we briefly reiterate why it is non-trivial 
to determine reliably the properties of phase transitions even in weakly coupled theories. 
Technically, this is due to the so-called infrared problem of thermal field theory [3]. Let g 
denote a renormalized gauge coupling. Then the theory has three different momentum scales: 

3 The conclusions of our analysis also apply to models beyond the Standard Model that resemble the 
considered effective theory at low energies, for instance certain triplet extensions of the MSSM [21] or the 
Standard Model with an extra coloured scalar |35| . 
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• The scale p ~ 2irT, technically originating from fields carrying non-zero Matsubara 
modes, is a "hard scale" ; perturbation theory at this scale is free from infrared problems, 
and is expected to converge well (practical tests in the Standard Model suggest that 
after a 2-loop computation the remaining errors are on the percent level [36 ). 

• The scale p ~ gT, called the Debye or "soft" scale and technically originating from the 
interactions between the hard modes and the long- wavelength excitations, is also free 
from infrared problems. In general, perturbation theory at this scale converges slower 
than that at the hard scale, but in the electroweak theory there are very many hard 
modes, meaning that the Debye scale comes with a large numerical coefficient, making 
it in practice not much smaller than the scale p ~ 2ttT. Indeed it has been checked (by 
comparing with full 4d simulations in an SU(2)+Higgs toy model [37]) that even the 
soft scale can be integrated out with good precision in the Standard Model [36] . 

• The "ultrasoft", or "magnetic" scale p ~ g 2 T/ir, first identified by Linde [3], is purely 
non-perturbative in nature, and needs to be studied numerically 

In the presence of a Higgs phase, the W^,Z a bosons have masses m w ,m z ~ gv(T)/2. If 
v(T) ~ T, than m w ,m z are formally "soft" scales and perturbation theory may or may not 
work, depending on numerical coefficients. However, if v(T) < gT, as is the case on the side 
of the high-temperature phase, then perturbation theory certainly breaks down. Therefore, 
in the following, m w and m z are treated as ultrasoft scales. 

Now, even though a scale hierarchy of the type described is a problem for perturbation 
theory (as it happens, it is also a challenge for direct 4d lattice simulations), it is a blessing 
once the problem is rephrased in an effective field theory language. Indeed the hard and soft 
scales can be integrated out perturbatively; it is only the ultrasoft scale which needs to be 
studied with lattice simulations. The integration out is called dimensional reduction [38| [39] . 
and the lattice simulations are then those of purely bosonic 3d gauge+Higgs theories. 

It is worth stressing that even if no simulations were carried out, it would nevertheless be 
useful to organize the computation in the above language. The reason is that the integrations 
over the hard and soft scales implement all-orders resummations (such as the "daisy" one), 
which are in fact necessary even for observables not sensitive to the ultrasoft scale. 

With this background, the general steps of the adopted approach are as follows: 

• Derivation of a 3d effective theory. The first step is to integrate out the hard 
and soft scales. To be specific, this step involves not only the actual finite-temperature 
calculations but also the corresponding vacuum computations, in order to express the 
MS scheme parameters in terms of experimentally measurable quantities. 

• 2-loop perturbation theory within 3d effective theory. Once an effective the- 
ory is available, it can first be studied with 2-loop perturbation theory. This is the 
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level that previous experience from the Standard Model and MSSM has shown to be 
semi-quantitatively accurate, provided that the transition is strong enough. Analytic 
expressions also yield a qualitative understanding of various parametric dependences. 

Lattice formulation of 3d effective theory. The derivation of the effective theory 
and its 2-loop perturbative analysis are, as a rule, carried out in a continuum regulariza- 
tion scheme, such as the MS scheme. Obviously, a lattice provides for a regularization 
scheme of its own. For a systematic study, the two schemes need to be related to each 
other; the principal relations have been worked out for a large class of theories |40j . 



• Numerical simulations within 3d effective theory. The last step of the program 
is to carry out lattice simulations. Although these are substantially simpler than full 
4d simulations, they do remain non-trivial: infinite-volume and continuum limits need 
to be carefully taken in order to obtain physical results. 



Without going into further details, which have been explained in refs. [261141]. we recall the 
continuum form of the 3d theory simulated. The theory contains two Higgs SU(2) doublets, 
Hi and H2, and a field U which is SU(2) singlet but SU(3) triplet. The Lagrangian has then 
the most general form allowed by symmetries, 

C 3d T* = ±TrC% j + (D!U?(D!U) + rr£(T*)tfU + \ ir (rfU) 2 

+ 7l U^UHlH-i + 7 2 U^U H J 2 H 2 + [ 7l2 U^U h\h 2 + H.c. ] 

+ + {DfH^iDfH,) + {DfH^{DfH 2 ) 

+ m 2 (T*) h\h x + m 2 (T*) h\h 2 + [m\ 2 {T*) H\H 2 + H.c. ] 

+ X^HlH-i) 2 + X 2 {H ] 2 H 2 ) 2 + Agfl'J^fl'Jfla + Xa h i H 2 H I h i 

+ [\{H\H 2 ) 2 + \ & H\H x h\h 2 + \ 7 hIh 2 h\h 2 + H.c. ] , (2.3) 

where a factor T* has been inserted in order to keep 4d dimensionalities of fields and couplings; 
Df, Df are the SU(3) and SU(2) covariant derivatives; G^, the corresponding field 
strength tensors; and H 2 = ia 2 H 2 . The hypercharge coupling has been set to zero, so there 
is only a global U(l) symmetry. The gauge couplings are denoted by g\, and g 2 for SU(2) 
and SU(3), respectively. The actual values corresponding to the setting of eq. (j2.2j) are 

m\(T*) « 26504 GeV 2 + 0.1311(T*) 2 , (2.4) 

m 2 2 (T*) w -4004 GeV 2 + 0.6311(T*) 2 , (2.5) 

m 2 12 (T*) fa -1481 GeV 2 - 0.0133(T*) 2 , (2.6) 

ml (T*) rj -4958 GeV 2 + 0.8607(T*) 2 , (2.7) 

7l « -0.0031 , 7 2 « 0.9995 , 7l2 « -0.0018 , X v w 0.2020 , (2.8) 
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volumes 






8 


12 3 , 


16 3 






10 


16 3 








12 


16 3 , 


20 3 


32 3 , 12 2 x 36, 20 2 


x 40 


14 


24 3 , 


14 2 


x 42, 24 2 x 48 




16 


24 3 , 


16 2 


x 48, 20 2 x 60, 24 2 


x 72 


20 


32 3 , 


20 2 


x 60, 26 2 x 72, 32 2 


x 64 


24 


24 3 , 


32 3 


48 3 , 24 2 x 78, 30 2 


x 72 


30 


48 3 









Table 1: The simulation volumes used in multicanonical simulations at temperatures around T* . 

A 2 ps 0.1491 , A 3 « 0.0589 , A 4 ps -0.1784 , (2.9) 
A 6 ps -0.00113 , A 7 « -0.00117 , (2.10) 
g 2 s « 1.085 . (2.11) 

As can be seen from eqs. (|2.5p . (|2.7p . it is the fields H 2 and U that "drive" the transition; 
correspondingly, the couplings 7 2 , A^, and A 2 are the most significant ones. 



Ai « 0.0649 , 
A 5 ps 0.00009 , 
Ol ~ 0.418 , 



3. Lattice simulations 
3.1. Action and algorithms 

The theory in eq. (|2.3p is discretized in the standard way, and for implementation details, 
we refer to refs. |26^ I41j . For the gauge action we use the usual single plaquette Wilson 
formulation, with SU(2) and SU(3) lattice couplings j3 w = 4/(g 2 y T*a) and (3 S = Q/(g 2 T*a). 
Here a is the lattice spacing, which we parameterize through j3 w from now on. 

Only the bare mass terms require renormalization in a 3d super-renormalizable theory: 

m^u = m 2 + Am 2 . (3- 1 ) 

Here m 2 stands for either m 2 , m 2 , m 2 or m 2 2 in the MS scheme, and Am 2 is the counterterm 
containing the linear and logarithmic in a divergences |26j . With this renormalization only 
0(a) cutoff errors remain, and the continuum limit can be taken by performing simulations 
at different /3 W and by extrapolating /3 W — > oo afterwards. 

The update algorithm is a combination of heat bath and overrelaxation updates, with one 
compound update sweep consisting of one heat bath and three overrelaxation sweeps. The 
measurements are performed and recorded every two compound update steps. 

It turns out that at the critical temperature the transition is strongly of first order, and 
the system does not spontaneously tunnel from one metastable phase to another with stan- 
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Figure 2: Left panel: the expectation values of (h\H 2 ) and (U^U) as functions of the temperature 
T* in the continuum limit. The vertical dashed line shows the phase transition temperature. Right 
panel: continuum extrapolation of (H\H 2 ) at four chosen temperatures. 



dard simulation algorithms. However, the probability distribution at the tunnelling region 
needs to be accurately resolved in order to precisely determine the critical temperature, the 
order parameter discontinuities, and the surface tension. We handle the strong metastability 
through the use of the multicanonical simulation method with automatic weight function 
calculation, combined with reweighting of the simulation temperature. We use the average 
h\h 2 as the multicanonical order parameter since this condensate is the most sensitive to 
the phase transition {H\ is practically inert due to the large value tan /3* = 15). The volumes 
simulated are listed in table [TJ 



3.2. Observables and results 

The observables measured on the lattice are all extracted from gauge-invariant operators and 
are therefore gauge-independent by construction. Here we reiterate their definitions and show 
the main results. 



Condensates as functions of the temperature: In order to obtain an overall view of 
the behaviour of the condensates {h\h^ and (U^U) we perform a series of simulations at 
T* = 70-90 GeV at three different lattice spacings, (3 W = 12, 16 and 24. The values of these 
condensates in the MS scheme with the scale parameter p, = T* are obtained by subtracting 
lattice divergences, for instance 



mH 2 ) 

(T*) 2 



H\H 2 ) 



latt 



(T 



g 3g 2 
2vraT* 16vr 2 



log 



6 



aT* 



+ 0.66796 



+ 0(a) , 



(3.2) 
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Figure 3: The volume- averaged probability distribution of H 2 ^2 a ^ three different volumes with 
fi w = 16 and at T*^ _ 16 . As the volume increases, the probability density between the two peaks 
decreases exponentially. 

where £ = 3.1759... originates from a 3d lattice tadpole integral. Here and in the following, 
we refer to the condensates with their original 4d dimensionalities. For our parameter settings, 
the values of these condensates extrapolated to the continuum are shown in the left panel 
of fig. [2J The figure highlights a first-order phase transition at a temperature around 80 
GeV (the transition temperature as determined with separate multicanonical simulations is 
indicated with the vertical dashed line; see discussion below). The stop field U is also affected 
by the transition, due to its strong coupling to B.2- 

On the right panel we show the continuum extrapolation at four selected temperatures. 
The lattice divergences have been subtracted according to eq. (|3.2p . whereas the remnant 
0{a) effects have been eliminated by a linear extrapolation in 1/ j3 w . 

Critical temperature (T*): The critical temperature is defined by the value of T* at 
which two phases, identified through the expectation value of (H2H2), are equally likely to 
exist. Because the tunneling between the two phases is strongly suppressed (it is exponentially 
suppressed at large volumes, cf. eq. (13.8P ). multicanonical simulations are implemented to 
overcome the tunneling barrier. 

In our case the critical temperature is expected to be at around 80 GeV, as already indicated 
by fig. [2j Near this temperature we therefore run multicanonical simulations for the volumes 
and lattice spacings listed in table[TJ For most of the lattice spacings, several volumes are used, 
enabling us to crosscheck the absence of finite-volume effects. The number of measurements 
at each volume is around (0.5 - 2) x 10 6 . Cylindrical volumes are needed for the surface 
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Figure 4: The critical temperature as a function of the lattice spacing 1/(3 W . With the volumes shown 
in table [T] no systematic volume dependence is seen, and all volumes are included in the plot. The 
interpolating curve is a second-order polynomial in a, fit to the points /3 W > 14. 



tension measurement, as described below. 

As an example, the distribution of the volume-average of H 2 H 2 , resulting from the simu- 
lations performed at f3 w = 16, is shown in fig. [3l Clearly, the probability of configurations 
between the two phase peaks becomes strongly suppressed as the volume increases. The 
temperature has been reweighted from the simulation temperature (T* = 79.5 GeV) to the 
apparent critical temperature where the area of the peaks is equal. The critical temperature 
averaged over all three volumes is T*„ _ 16 ~ 79.57 GeV. 

All the critical temperatures obtained for the lattice spacings and volumes of table Q] are 
summarized in fig. SJ A clear dependence on the lattice spacing is present, but results tend 
to stabilize as the continuum limit is approached. We then fit only data at (3 W > 14 by a 
second-order polynomial in 1/(5 W . The continuum intercept reads 

T* = 79.17 ± 0.10 GeV , (3.3) 

with x 2 /d-o.f. = 11.2/8 for the fit. 

In fig. [5]the joint probability distribution of the volume-averaged H 2 H 2 and U^U is shown. 
Clearly, the condensates move together: when H 2 H 2 becomes large, this effectively increases 
the mass of the [/-field through the interaction mediated by the coupling j 2 111 e Q- (|2.3j) , and 
therefore its fluctuations become smaller. In fact in the MS scheme with the scale choice 
p, = T* the value of (U^U) becomes slightly negative at low temperatures, but this simply 
means that U is tightly confined into its "symmetric" phase in this regimeu 

4 It is in principle possible that separate symmetric f-> broken H2 and U phase transitions exist, giving rise 
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Figure 5: Contour plot of the joint probability distribution of the volume-averag ed HlH 2 and WU 
at T*. The measurement is for (3 W = 20, volume = 32 2 x 64. The MS scheme value of (WlJ) is slightly 
negative in the low-temperature phase (cf. the text). 



Higgs discontinuity (v(T*)): The Higgs discontinuity is defined on the lattice in a gauge- 
invariant and scale-independent manner as 

*M) =A(£aJfl,). (3-4) 

/ latt i=l 

The quantity A(-) = (•) broken — (^symmetric is measured from the probability distributions 
at the critical temperature (such as those shown in fig. [3]) by integrating over the peaks 
independently^] In our case, the contribution of H\ can be neglected since tan/3* S> 1. 

The Higgs discontinuity as a function of the lattice spacing is shown in fig. [6] (left panel) . 
Each point is obtained by averaging over the volumes listed in table [U In this case v(T*)/T* 
levels off almost completely at small enough lattice spacing, and a linear fit to points with 
Pw > 16 gives 

v(T*) 

- 1.117 ±0.005 , (3.5) 



T. 



c 



with x 2 / d -a f - = 2.8/3. 



to three phases: fully symmetric, broken Hz, and broken U. This has indeed been observed in a different 
parameter region, see fig. 9 of ref. [26]. However, in the case at hand we do not observe a stable broken U 
phase (no scan in has been performed here). 

5 This involves setting a "separatrix" between the two phases, which is chosen to lie at the minimum of 
the probability distribution. As can be deduced from fig. [3] the ambiguity related to the choice becomes 
exponentially insignificant at large volumes. 
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Figure 6: The Higgs discontinuity (left panel) and the latent heat (right panel) as functions of the 
lattice spacing. The linear continuum fit to points at (3 W > 16 (solid line) is also shown. 



The value of v(T*)/T* is large enough so that the sphaleron rate (which is proportional 
to the baryon number violation rate) after the phase transition is negligible. Nevertheless, 
it should be stressed that v(T*) as denned by eq. (|3.4p is a purely thermodynamic quantity, 
and therefore is not equivalent in any strict sense to the sphaleron rate. However, as non- 
perturbative real-time simulations for the Standard Model have shown [32], the value of 
v(T*)/T* is strongly correlated with the sphaleron rate. Given the fact that non-perturbative 
real-time simulations are very expensive, only v(T*)/T* is measured in the present study. 



Latent heat (L) : The latent heat is defined as the discontinuity of the energy density across 
a first order phase transition. It plays an essential role in the real-time hydrodynamics of 
bubble nucleation and growth. Within the approximate parametrization of the 3d theory used 
in this work, in which only mass parameters depend "non-conformally" on the temperature 
[cf. eqs. ([2.4p - (|2.1ip ], it can be measured as 



*\3 



\ dT* 

+ (h\h. 



m 2 (T* 



(T 



*\2 



dT* 



(T*) 2 



d 

dT* 



m? 2 (T*) 



H.c. 



(t*) 2 j 1 ~~7 / • (3-6) 

The discontinuity in eq. (|3.6p is readily measured from multicanonical simulations, with the 
results shown on the right panel of fig. [6j The continuum limit gives 

L 



*\4 



0.443 ± 0.004 



(3.7) 



with x 2 /d.o.f = 1.01/2 using a linear fit to data with f3 w > 16. 
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Figure 7: The surface tension of the interface between the symmetric and broken Higgs phases, 
plotted versus the inverse area of the interface. 



Surface tension (<r): The surface tension is defined as the additional free energy per 
area carried by an interface between the two co-existing phases. This means that, in the 
large- volume limit, the probability of a configuration which contains an interface of area A, 
denoted by P m i n , is smaller than the probability of a configuration without interfaces, -P m ax> 
by a factor 

* «P(-£V (3-8) 



P V T 1 * I 

■'max \ J. / 

When the volume is finite, there are corrections to this relation, and properly accounting for 
these accelerates the convergence to the infinite- volume limit (cf. refs. \26\ W\\ for details). 

The probabilities P min and P max can be directly read from the distributions in fig. [3j P max is 
the average peak height and P min is the minimum between the peaks. Because of the periodic 
boundary conditions the configurations here contain (at least) two interfaces. In practice it 
is advantageous to use lattice volumes where one dimension is longer than the other two, e.g. 
cylindrical volumes (L z ^> L x , L y ). In this case the interfaces are oriented transversely to the 
long direction (A sa 2L x L y ). 

In fig. [7] measurements of the surface tension at each of the cylindrical volumes and lattice 
spacings indicated in table Q] are shown. In this case we are not able to obtain independently 
reliable (i.e. with small x 2 /d.o.i) infinite volume and continuum limit extrapolations. This is 
likely due to still remaining finite volume effects; surface tension measurements are notoriously 
sensitive to volume. However, it is clear from the plot that the measurements settle down on 
a narrow band, independently of the lattice spacing. Extrapolating the band to the infinite 
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m* = 7 TeV, m!,= 70.5 GeV, nC = 126 GeV 

Q U h 




Figure 8: Comparison of pcrturbativc and lattice results for the properties of the phase transition 
(here v(T*) refers, strictly speaking, to different quantities on the two sides; cf. sec. 14. 2[) . 



volume we cite a conservative but unprecise estimate of the error: 

a 



-= 0.035 ± 0.005 . (3.9) 



4. Comparison with perturbation theory 

With a view of learning about generic features of the dynamics of the theory, probably appli- 
cable also to other parameter values than the very ones considered here but nevertheless close 
to rrih ~ 126 GeV, we proceed to comparing the lattice results with those of 2-loop pertur- 
bation theory within the 3d theory. We stress that since both results are based on the same 
3d theory, the comparison is not jeopardized by perturbative uncertainties in dimensional 
reduction and vacuum renormalization as discussed in sec. 12.11 Indeed, these ultraviolet fea- 
tures play a role only in the relation of the approximate parameters (T*, m~ , etc) to the 
physical ones (T, m/j,m^ fl , etc). For conceptual clarity, we furthermore split the comparison 
into two parts, given that some of the perturbative numbers cited are specific to Landau 
gauge, in accordance with established (although not necessary) conventions of the field. 

4.1. Identical observables 

Two of the observables, namely the critical temperature and latent heat, have definitions 
[see eq. (|3,6p for the latter] that can be operatively applied both to lattice and perturbative 
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calculations. Their values for the case analyzed are 

79.17(10) (latt) , = 84.4 (pert) , (4.1) 

0.443(4) (latt) , = 0.26 (pert) . (4.2) 

For the critical temperature the situation is illustrated in fig. [8l As has been observed also in 
the past [26], the main qualitative effect from non-perturbative dynamics is that the critical 
temperature is lowered. The latent heat is enhanced by ~ 50%. 

4.2. Correlated observables 

Within 2-loop perturbation theory, the gauge-independent observable defined by eq. (|3,4p 
happens to be very close to the gauge-fixed Higgs vev as computed in Landau gauge. Due to 
the fact that the Landau-gauge convention continues to be widespread in the literature, we 
therefore compare the lattice number directly with the Landau-gauge perturbative result: 

= 1.117(5) (latt) , = 0.9 (pert) . (4.3) 

eq. I|3.4|) \ J c / Landau 

The percentual strengthening effect is smaller than for L, because L is essentially quadratic 
in v (cf. eq. (|3.6p ). Another quantity for which we are influenced by convention and ease of 
computation is the surface tension; this is usually extracted from the Landau gauge effective 
potential, with tree-level kinetic terms employed in finding the saddle point solution (for a 
recent discussion, see ref. [13]). The comparison reads 

(7^3) = °- 035 ( 5 ) ( latt ) > ( 7t£w) = °- 025 (P ert ) • ( 44 ) 

\\ 1 c) / eq. Vl J cJ / Landau 

5. Discussion and conclusions 

The recent LHC discovery of a Higgs- like boson with a mass of around 126 GeV may have pro- 
vided crucial information for electroweak baryogenesis. In many models beyond the Standard 
Model, the success of electroweak baryogenesis in explaining the baryon asymmetry of the 
Universe is indeed very sensitive to the Higgs mass through the requirement of a strong first- 
order electroweak phase transition. In the MSSM, perturbative studies have suggested that 
a strong first-order electroweak phase transition may exist even at ~ 126 GeV \22\ 125] . 
However, since the transition is not exceedingly strong and the side of the "symmetric" phase 
is purely non-perturbative in nature, it is not clear whether the perturbative predictions are 
quantitatively accurate. In this paper we have studied the infrared dynamics of the transi- 
tion by simulating a dimensionally reduced effective theory numerically, and compared the 
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L 
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results with 2-loop perturbative calculations within the same effective theory. Unless there 
are unexpectedly large 2-loop corrections to the relations between the effective parameters 
of the dimensionally reduced theory and four-dimensional physical low-energy observables, 
the simulations correspond to an MSSM-like parameter point with rrih ~ 126 GeV and 
m~ ~ 155 GeV, with larger uncertainties on the latter. (The simulations, however, are 
expected to also cover other extensions of the Standard Model having a stop-like field at the 
electroweak scale, i.e. providing for a similar low-energy effective theory.) 

In the lattice simulations carried out, we have consistently seen a stronger transition than 
in perturbation theory. Actually, despite the larger Higgs mass, the strengthening effect is 
more substantial than in ref. [26]. In some sense the system is driven towards a phase where 
the right-handed stop experiences very strong fluctuations (manifested by a large (U'U)), 



and the transition to the electroweak minimum takes place "from there" (cf . fig. [5j) |£| 

For a precise understanding of baryogenesis, it is not enough to study the properties of 
the transition at the critical temperature, but issues such as supercooling, nucleation, and 
bubble dynamics need to be considered as well (see e.g. refs. [45j-|47| for recent discussions). 
The nucleation temperature, T n , is well approximated by the classical estimate [48] if it is 
calculated with non-perturbative values of the latent heat and surface tension inserted [49] , 
After nucleation and bubble collisions, the latent heat released may also reheat the system 
towards the critical temperature, which would enhance baryon number washout. 

Various scenarios for the real-time dynamics of the transition have been studied in ref. |50j . 
and in fact the dynamics of the present transition is not unlike case (A) considered there. 
More precisely, reheating up to T c would take place if 



This is (narrowly) avoided according to eqs. (|4.2p and (|4.4p . Moreover, supercooling is roughly 



It is therefore quite modest, and is not expected to change the Higgs vev substantially (cf. 
figs. O [8]). However, modern hydrodynamics studies of the phase transition will make these 
conclusions firmer |51j . 

At the moment the physical 4d parameter values to which our simulations correspond, 
particularly the right-handed stop mass, contain uncertainties of several GeV. To remove 
this perturbative uncertainty, full 2-loop dimensional reduction and at least 1-loop on-shell 
vacuum renormalization computations (expressing MS scheme parameters in terms of physical 
low-energy observables), such as were carried out for the Standard Model |33j . are needed. 

6 Considerations such as those in ref. [44] are evaded because this is not a perturbatively "colour-broken" 
phase and the transition to the physical vacuum does take place, as the estimate in eq. Q5.2p shows. 
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Despite the encouraging results that we have found, it is also clear that an exclusion of 
light SU(2)-singlet stop-like particles at the LHC could easily rule out the MSSM-based 
electroweak baryogenesis scenario. Therefore, it may be worthwhile to apply the techniques 
of the present study to more general models. Very many possibilities can be envisaged; as 
an example of a relatively well-constrained one, let us mention the so-called Inert Doublet 
Model \52\ 153] . i.e. a particular version of the two-Higgs-Doublet Model with an imposed 
unbroken Z(2) symmetry which reduces the number of free parameters. This model has 
many attractive features, for instance the heavy Higgs doublet could naturally serve as Dark 
Matter [541 [55] . The theory could conceivably lead to a strong first order phase transition as 
well [56l EZl [58] , and a non-perturbative study may again be welcome. 
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